Correlation functions for Id interacting fermions with spin-orbit coupling 
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We compute correlation functions for one-dimensional electron systems which spin and charge de- 
grees of freedom are coupled through spin-orbit coupling. Charge density waves, spin density waves, 
singlet- triplet- superconducting fluctuations are studied. We show that the spin-orbit interaction 
modify the exponents and the phase diagram of the system, changing the dominant fluctuations and 
making new susceptibilities diverge for low temperature. 
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In past two decades there have been an intense ef- 
fort in studying quasi-one-dimensional electron systems 
(Q1DES). This interest has its origin in the simplicity 
of the models which describe them, and at the same 
time, in the possibility of making contact with experi- 
ments. Examples of these Q1DES of recent construction 
are carbon nanotubesi, conducting polymers^ and semi- 
conductor heterostructures&. From the theoretical point 
of view the simplest formulation of a Q1DES is given 
by the Tomonaga-Luttinger^ 5 . model which describes the 
major qualitative features of the interacting Q1DES such 
as the spin-charge separation and the non universal ex- 
ponents in the decay law of correlation functions. 

In realistic situations the electrons are moving in elec- 
tric fields inside the materials: microscopic and macro- 
scopic ones, the latter responsible of confining the elec- 
trons to a reduced region of space. As a consequence it 
appears a magnetic field in the rest frame of the electron 
which couples with its intrinsic magnetic moment and 
breaks the spin-rotation SU(2) symmetry. This is known 
as spin-orbit (SO) interaction (or spin-orbit coupling). 
Despite its relativistic origin, this interaction has an im- 
portant effect on existing two-dimensional electron gases 
(2DEG) such as GaAs/AlGaAs^ 7 . and InGaAs/InAlAsM 
heterostructures. It is responsible for the modification 
of their band structure by lifting the spin degeneracy 
and for positive magnetoresistance effects, known as weak 
antilocalizationiS 

In Q1DES there exists an additional potential, respon- 
sible of confining the electrons in a narrow channel pat- 
terned in 2DEG heterostructuresii. Although as far as 
we know there is no experimental evidence or measures 
of the strength of the SO coupling resulting from such a 
potential, theoretical work indicates that it affects quan- 
titatively the splitting energy behavior as a function of 
the wave vector k, and produces an asymmetric defor- 
mation of each spin branch, i.e. the Fermi velocities take 
different values for different directions of motion^. 

In Ref. |13| the following Hamiltonian was proposed as 
a model for a Q1DES with SO coupling 



H — H + Hi nt 

where the noninteracting Hamiltonian is 



(1) 



H = vi [(& + ki)c\ m c km - (k - k 1 )cl Ll c kLl 



V2j2[(k + k 2 )cl Rl c kRl -(k-k 2 )cl L1 c kL1 . (2) 
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It consists of a modified Tomonaga-Luttinger model that 
takes into account the asymmetry in the spectrum for 
each spin branch making vi ^ v 2 (and k\ ^ k 2 ). cjL cre- 
ates a right-going (r=+l) or left-going (r=-l) electron. 
The interacting Hamiltonian describes forward scattering 
electron-electron interactions and has a standard form 5 . 
Umklapp and backscattering terms are irrelevant if wc 
are far from half Ofilling in the former case, and we re- 
strict to repulsive interactions in the latter. 

In this article we compute correlation functions for 
charge density wave (CDW), spin density wave (SDW) 
4fcp charge density wave (4fcp) and singlet- and triplet- 
superconductivity (SS and TS) operators for the model 
presented above. The correlation functions for these 
operators are well-known in the case of zero spin- 
orbit coupling 5 ^^, including logarithmic correction 
factor o 16 ! 17 and time and temperature dependence^ 5 -. We 
extend these calculations to the case in which SO inter- 
action are present, and study how the exponents of their 
algebraic decay are modified. We find interesting mod- 
ifications of the phase diagram of the system when SO 
interactions are present. For certain regions of the pa- 
rameter space, SO coupling changes the dominant fluc- 
tuations, and makes new susceptibilities diverge for low 
temperature. 

The Hamiltonian ffl can be studied by the use of 
bosonization techniqu o 5 i 18 as in Ref. 0. For convenience 
we shall define an average velocity vq = (vi + v 2 )/2 and 
the difference Sv = v 2 — v% , and the same for the Fermi 
momentum fco = (&i + k 2 )/2 and 5k = k 2 — k\. If wc 
introduce the usual phase fields 4> p ((f> a ) for charge (spin) 
degrees of freedom and the dual field n p (n CT ) the Hamil- 
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tonian can be represented as 



H =-£ I dx 



— (d^ p ) 2 +K p U 2 p 



+ - 



dx 



+Sv / dx [{d x( /) p )u a + (d x ^)n p ] 



(3) 



are the propagation velocities of the spin and charge 
collective modes of the decoupled model (Sv = 0), and 
K pa are the stiffness constants. The spin-orbit inter- 
action appears as an effect that breaks the spin charge 
separation as reveals the presence of the third term in the 
last equation. Nevertheless the Hamiltonian © can still 
be diagonalized in terms of two new phase fields which 
contains a mixture of spin and charge degrees of freedom. 
The propagation velocities of these collective modes are 



J p ,a 



transition^. For K p > K a , 5v p turns out to be lower 
than Sv a and as Sv grows from the zero value the physical 
divergence takes place in the charge compressibility. This 
instability is known as phase separation and has been 
shown to occur in the extended Hubbard models and in 
the t-J modelSi. In the case that K p < K a) the instabil- 
ity takes place in the spin subsystem and is related to the 
so called metamagnetic transition, observed for instance 
in the quasi-one-dimensional compound Ba 3 Cu 2 04C12 22 . 
It also arises in the phase diagram of the XXZ model 
with next-to-nearest neighbors^. In presence of a chem- 
ical potential (magnetic field), the region where k (x) is 
negative is associated with the coexistence of two phases 
with different hole concentration (magnetization). The 
divergence of k was found in other models with asym- 
metric dispersion 2 *. 

Let us now focus our attention on the correlation func- 
tions. Our interest in this work is to obtain their space- 
time and temperature T = 1//3 behavior. The operators 
for CDW, SDW, 4k F , SS and TS fluctuations in their 
bosonized form arc 



+ Sv 2 




Sv 2 



V 2 +V 2 + V p V a 
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As Sv —>■ 0, v + —+ max(u p , v a ) and w_ — * min(w p , v a ). As 
Sv increases, w_ decreases until it vanishes at the points 
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&4k F 

0SDW,z 
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— cos(2k,QX + V2ir(f) p ) cos(5kx + VZkOo- 



(9) 
(10) 



(11) 
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(5) 
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cos(2fcoa; + v 2-K<j) p ) sm(5kx + v2-n:6 a ) 



At these points, the freezing of the lower bosonic mode 
is accompanied by a divergence in the charge and spin 
response functions. The static charge compressibility k 
diverges at Sv — Sv p , and at Sv = 5v a occurs a divergence 
of the static spin susceptibility x- They behave as 
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(13) 
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where kq and xo are the values of k and x m absence 
of SO coupling. Beyond these points the susceptibil- 
ities becomes negatives. This behavior of the static 
response functions together with the vanishing of the 
collective modes velocity indicates that the system be- 
comes unstabl e^ 1 ? and undergoes a first order phase 



where a is a short distance cutoff and 6 V is related to the 
conjugated field by the relation II \ = d x 9\. 
The correlation functions are defined as 



Ri(x,T;(3) = (t t 1 (x,t)01(Q,O) 



(17) 



where T T is the (imaginary) time-ordering operator. 
These objects were calculated in the path integral frame- 
work within the Matsubara imaginary time formalism 
and the results are: 
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Rcdw{x, r; (3) = Rsnw,z(%, t; (3) = 

cos 2k$x 



2(7ra) 5 



+K a vZ)/2 



Z + Z- 



H sign(a;r ) 



h.c. 



(18) 
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COsAk X _ 2K p vl , - ,-2K p 

R 4kF (x,T;l3) = — — — (z+z+) p +(z_z_) 



(20) 



R$${x, t; /3) = Rxs >0 (x, r; /3) = — — — (z + z + ) 



(21) 



i?TS, ±1 (x,r;/3) = ^(z + z- + )-K^ + ^/^)/ 2 (z_z_)-^^ + ^^)/ 
(27ra) z 



±G signer) 



(22) 



where 



sin ^( w ±l r l + e + "0 
sin-^% 



sin 



(23) 
(24) 



and the exponents depend on the the stiffness constants 
multiplied by the factors that include mode velocities de- 
pendence. They are given by 



*4=± 



v± 



9 2 

v+ — v_ 



(25) 
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5v u± - (<5w5J - Sv 2 ) 



v± v\ — v_ 



with A = p, a, and 



H =Sv 



G=Sv 



KpVp + K a v a 



Vp/K p + VajKo 



(27) 

(28) 
(29) 



v± and [x± are positive, and 8±, G and H has the same 
sign as Sv. 



In the model with zero SO coupling the SU(2) symme- 
try can be restored by imposing the constraint K a = 1, 
which emerges naturally if the model under study is the 
continuum limit of a lattice model with only charge den- 
sity interactions. In this case this is not possible; the 
SU(2) symmetry keeps broken even for K a = 1 as reveals 
the differences in the decays between SDW operators cor- 
relation functions in the z direction and in the x, y direc- 
tions. As in the zero SO case, correlation functions for 
SDW operators in the z direction and CDW operators are 
equal, and the same happen with TS,0 and SS operators. 
This degeneracy is broken by logarithmic corrections that 
arises if irrelevant backscattering or umklapp terms are 
included^. 

An interesting point to observe is the appearance of 
two terms in the SDW, xy correlation functions (Eq. I19|) 
where the modulations have different frequencies and de- 
cay with different exponents. As 6± has the same sign 
as Sv (see Eq. [I7\ and the comment below Eq. I29fl for 
V2 > v\ (v% < vi) the dominant term is the one with 
frequency &2 In other words the biggest frequency 

dominates. Also i?TS,±i becomes oscillating. 

Up to here we have obtained very general formulae for 
space-time and temperature dependent correlation func- 
tions for the model under analysis. We can gain physical 
insight by observing the algebraic decay of the instan- 
taneous correlation functions at zero temperature, and 
studying how the exponents get modified from the zero 



4 



a 





CDW SDW SS 


TS 




1 














„ - Q 72 


0.4 0.6 


N . i 


1 






\ 


2 






\ 

\ 



5v 



FIG. 1: Behavior of the exponents a' t s as a function of Sv 
(in units of vo)- For n p = 1.2vo, v„ = 0.8uo, K p = 0.6 and 
T^ct = 0.85. For Sv > 0.16 SDW, xy fluctuations get dominant, 
and for Sv > 0.25 ass gets positive, and xss divergent for 
T -> 0. 



SO case. The functions behave as 



Kp 




FIG. 2: Phase diagram in K p — K a space. The phase in 
brackets is the subdominant one, which becomes dominant 
for strong enough SO coupling. 



Ri(x) 



\-2+ai 



(30) 



The exponents a^s determine the divergence of the corre- 
sponding Fourier space susceptibility as T — > 0, Xi{T) ~ 
T~ a &. This makes these instabilities of a completely 
different nature than the ones described in equations l(7)l 
and JHJ. The expressions obtained for the cti are 
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"cdw = "sdw,z =2 - K p v p - K a v a (31) 

asnw, x = a S DW, y =2(1 + \9°\) - K p v p - ^/K a (32) 

ass = «ts,o =2(1 + \0 p \) - ii p IK p - K a v a (33) 

axs,±i =2 - tiP/Kp - v?lK a . (34) 



These are the new exponents, which retain the same 
structure as in the zero SO coupling, but modified by 
the factors 
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/i A = /4 + ^ (35) 
v x = v\ + v x _ (36) 
6 x = el + 9 x _. (37) 

When Sv — * 0, 8 X — > and fi x ,v x — > 1, so we reproduce 
the right results for the zero SO case. 

For finite SO coupling, Sv appears as a parameter 
which plays a role in determining the slowest decaying 
correlation function, and which are the divergent suscep- 
tibilities. In Fig. ^we observe as an example, the behav- 
ior of the exponents as a function of Sv for v p = 1.2t;o, 
v a = 0.8v , K p = 0.6 and K a — 0.85. For Sv small, 
CDW fluctuations are dominant, but for Sv > 0.16t>o 
the SDW,xy correlations decay slower. For small Sv, 
CDW and SDW are the only diverging susceptibilities 
for T —> 0, but for Sv > 0.25vq, ass becomes positive, 



FIG. 3: Phase diagram in K p — Sv space for v p = 1.2vo, 
v a = 0.8«o and different values of K a - Sv > 5v a below the 
doted line and metamagnetism occurs. 



and xss divergent for T —* 0. Calculations of the elec- 
tron band structure modified by SO coupling show that 
these values of Sv should correspond to typical Q1DES 12 . 

A careful analysis of the exponents allows us to con- 
struct a phase diagram in K p —K a space (Fig. 0). In each 
region we indicate the dominant fluctuation for small Sv, 
and in brackets, the dominant one for stronger Sv. Other 
subdominant fluctuations are not indicated. Cross sec- 
tions of the phase diagram are shown in Fig. [3] In this 
plot the K p — Sv space can be observed for K p < 1 and 
different values of K a . For small Sv CDW fluctuations 
dominate and for stronger Sv the system can be either 



5 



in the SDW or in the SS phase depending on the val- 
ues of K p and K a . In the region below the doted line, 
Sv > 5v a , the static spin susceptibility becomes negative 
and metamagnetism takes place 

In conclusion, we have computed correlation functions 
for a model of one-dimensional correlated electrons with 
SO coupling. This coupling destroys the spin charge sep- 
aration as was shown in Ref. Il3l and modifies the expo- 
nents of correlations decay. As a consequence the phase 
diagram gets modified. For strong enough SO coupling, 
it changes the dominant fluctuation, and makes new sus- 
ceptibilities diverge for T — ► 0. How logarithmic correc- 
tions originated in irrelevant backscattering and/or umk- 



lapp terms modify these results is an interesting problem, 
subject of future work. 
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